Preparation of arbitrary qutrit state based on biphotons 
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The novel experimental realization of three-state optical quantum systems is presented. We use 
the polarization state of biphotons, propagating in single frequency- and spatial mode, to generate an 

' arbitrary qutrits. In particular the specific sequence of states that are used in the extended version 

of BB84 quantum key distribution protocol was generated and tested. We experimentally verify the 
orthogonality of the 12 basic states and demonstrate the ability of switching between them. The 

' tomography procedure is applied to reconstruct the density matrices of generated states. 

o : 
fc: 
r> : 

(NJ . I. INTRODUCTION 

A. Three-state systems in Quantum Information 

<N '■ 

Q\ . The art of quantum state engineering, i.e., the ability to generate, transmit and measure quantum systems is of 
t-H ' great importance in the emerging field of quantum information technology. A vast majority of protocols relying on 
the properties of two-level quantum systems (qubits) were introduced and experimentally realized. But naturally, 
there arose a question of an extension of dimensionality of systems used as information carriers and the new features 
that this extension can offer. The simplest extension provokes the usage of three-state quantum systems (qutrits). 
Recently new quantum key distribution (QKD) protocols were proposed that dealt specifically with qutrits 0, Q 
and the eavesdropping analysis showed that this systems were more robust against specific classes of eavesdropping 
attacks 0, 3 • The other advantage of using multilevel systems is their possible implementation in the fundamental 
tests of quantum mechanics [jj , giving more divergence from classical theory. The usage of multilevel systems also 
provides a possibility to introduce very specific protocols, which cannot be implemented with the help of qubits such 
as Quantum Bit Commitment, for example [(| . Recent experiments on realization of qutrits rely on several issues. 
In one case, the interferometric procedure is used, where entangled qutrits are generated by sending an entangled 
r* ■ photon pair through a multi-armed interferometer [?J • The number of arms defines the dimensionality of the system. 
. £^ [ Other techniques rely on the properties of orbital angular momentum of single photons 0, 0, and on postselection 
of qutrits from four-photon states |lfj|| . Unfortunately all mentioned techniques provide only a partial control over 
a qutrit state. For example in a method, mentioned in 0, IS a specific hologram should be made for given qutrit 
state. The real parts of the amplitudes of a qutrit, generated in pj are fixed by a characteristics of a fiber fritter, 
making it hard to switch between the states. Besides, in this method no tomographic control over generated state 
had been yet performed. 



^1- 

o 



B. Biphoton properties 



In this paper we report the experimental realization of arbitrary qutrit states that exploits the polarization state 
of single-mode biphoton field. This field consists of pairs of correlated photons, is most easily obtained with the help 
of spontaneous parametric down-conversion (SPDC). By saying "single-mode" we mean that twin photons forming a 
biphoton have equal frequencies and propagate along the same direction. A pure polarization state of such field can 
be written as the following superposition of three basic states. 

|*) = ci|2,0) +c 2 |l,l) +c 3 |0,2), (1) 
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where c, = \ci\e % ^ i are the complex probability amplitudes. Bases states are the Fock states with certain number of 
photons in two orthogonal polarisation modes. It is reasonable to use vertical and horizontal polarisation modes for 
basis state definition. So by state |2,0) we mean that two photons have horizontal polarization, |0,2) is the basis 
state were both of the photons are polarized vertically and |1, 1) is the basis state with one horizontally polarized 
photon and one vertically polarized photon. Such a states can be generated in SPDC processes, first two via type I 
phase-matching SPDC and the last one via type II phase-matching SPDC. 

There exists an alternative representation of state |\&) that maps the state onto the surface of Poincare sphere [TT| 



l*> 



at(0»at(0',<//)|mc} 



and a'(6',(j)') are the creation operators of a photon in the certain polarization mode a^(6, < 



(2) 



at are photon creation operators in correspondingly horizontal and vertical po- 



where <r(6, < 

cos(9 /2)aJ. + exp(i(p) sin{6 /2)a^, aj., 
larization modes, 6,6' S [0, ir], 4>,(f>' S [0, 27r] are polar and azimuthal angles that define the position of each photon 
on a surface of a sphere. The values of the angles can be calculated using amplitudes and phases of Cj . 

The operational orthogonality criterion for the polarization states of single-mode biphotons was proposed in |l3j 
and experimentally verified in [14| . Let us consider setup shown on the Fig. ^ m order to formulate the criterion. 
Presented setup is a Brown- Twiss interferometer with a quarter-wavelength plate and an analyzer in each arm. Both 
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FIG. 1: Brown-Twiss interferometer with polarization selection. 



the plate and the analyzer in one arm could be rotated independently by an arbitrary angle. According to these 
angles there is a certain polarization of photons that corresponds to their lossless propagation through the arm. Let 
us assume that these polarization photon states are a' (6, (j))\vac) for transmitted arm and af(0', (/)') \vac) for reflected 
arm. In this case we would say that setup is "tuned" to detect the biphotons in polarization state |\&) 

If we have the biphotons in state in the input of the setup the coincidence count rate would be non-zero. The 
most interesting question is what happen with coincidence count rate R c when the biphotons in the input of the 
setup that is tuned to state |\&) are in any other state |^'). The answer was obtained in |L| where was shown that 
Rc ~ | (^1^') | 2 - So this scheme is a realization of the projector of an input state on the fixed one. In particular if 
we have input state orthogonal to the one setup tuned there would not be any coincidences in the scheme. Thus the 
coincidence rate is zero if and only if the input state is orthogonal to the one setup tuned. The significance of this 
criterion would be shown at the next section. 

The goal of our work was to demonstrate the ability to prepare any given polarization state |^) and as a straight- 
forward and practical example of given states, we chose the specific sequence that was presented in |l| . This sequence 
of 12 states can be used in an extended version of BB84 QKD protocol for qutrits. As it was shown one can construct 
4 mutually unbiased bases with 3 orthogonal states in each from any three orthogonal states using discrete Fourier 
transformation. Using a = |2, 0), (3 = |1, 1), 7 = |0, 2) as the first orthogonal basis we observed all 12 states shown in 
the table I. where a, f3 and 7 with fixed number of primes are the orthogonal states belong to the same basis and a 
scalar product of any two states belong to different bases is equal to 1/3. 



II. EXPERIMENT 



A. Experimental setup and method of measurement 



Our setup consists of two parts: first one performed the preparation of the biphoton field in a desired state and 
second one provide the measurement of the obtained states. 
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TABLE I: 12 states used in qutrit QKD protocol. 
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The preparation part of our setup (Fig. |2| is built on the base of a balanced Mach-Zehnder interferometer (MZI) 
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FIG. 2: Experimental setup (preparation part). 

The pump part consists of frequency doubled " Coherent Mira 900" femtosecond laser, operated at central wavelength 
of 800 nm, 75 MHz repetition rate and with a pulse width of 100 fs, average pump power was 20 mW. The Glan- 
Tompson prism (GP), transmitting the horizontally polarized fraction of the UV pump and reflecting the vertically 
polarized fraction, serves as an input mirror of MZI. The reflected part, after passing the compensation BBO crystal 
and a half- wave plate (HWP2), pumps two consecutive 1 mm thick type-I BBO crystals whose optical axis are oriented 
perpendicularly with respect to each other. The biphotons from these crystals pass through a 10 mm quartz plate 
(QPl) that serves as a compensator, and the pump is reflected by an UV mirror. Then the biphotons arrive at a 
dichroic mirror (DM) that is designed to transmit them and to reflect the horizontally polarized component of the 
pump coming from the upper arm of MZI. A piezoelectric translator (PZT) was used to change the phase shift of the 
horizontal component of the pump with respect to the one propagating in the lower arm. The UV beam, reflected 
from DM serves as a pump for 1 mm thick type-II BBO crystal. Two 1 mm quartz plates (QP2) can be rotated along 
the optical axis to introduce a phase shift between horizontally and vertically polarized type-I biphotons, and a set of 
four 1 mm thick quartz plates (QP3) serves to compensate the group velocity delay between orthogonally polarized 
photons during their propagation in type II BBO crystal. 

The measurement setup (Fig. |3J) consists of a Brown- Twiss scheme with a non-polarizing 50/50 beamsplitter; each 
arm contains consecutively placed quarter- and half waveplates and an analyzer that was set to transmit the vertical 
polarization. 

This sequence of waveplates and analyzer is referred to as a polarization filter as soon as it extracts a single-photon 
polarization state defined by their orientation. Interference filters of 5 nm bandwidth, centered at 800 nm and pinholes 
are used for spectral and spatial modal selection of biphotons. We use EGG-SPCM-AQR-15 single photon counting 
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modules as our detectors (Dl and D2). We should mention, that due to the low pump power, the stimulated processes 
in our setup are negligibly small and only pairs of photons have been generated. The measurement of the generated 
states is done using the tomography protocol that was developed for polarization qutrits p"?l Il8| . From experimental 
point of view we have to measure all fourth-order correlation moments in order to reconstruct coherence matrix 



K A = 




which contains the following moments 
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B= (&t b 2 ), C=(aWab), 
D = {afab), E = (a^b 2 ), F^{a%^b 2 ). 
In the experiment one can measure directly only the diagonal elements of the coherency matrix. For example one 
measures B moment when there are only two vertical polarizers in the arms and no plates , adding quarter wavelength 
plate rotated at 45° to one of the arms - C moment, adding quarter wavelength plates rotated at 45° to each arm - 
A moment. In all other cases the coincidence rate would be proportional to linear combination of different moments. 
By choosing an appropriate combination containing as few number of moments as possible was proposed the following 
tomography protocol shown on the table 2. 



TABLE II: Tomography protocol. Each line contains orientation of the half (9 s ,i) and quarter (Xs,i) wave plates in the 
measurement part of the experimental setup. Last column show the corresponding moment. 
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As a result of the series of 9 measurements with different orientation of the plates in the arms of Brown- Twiss 
interferometer we can reconstruct the coherency matrix. Since we deal with a three-dimensional Hilbert space density 
matrix of measured state determined by 8 real numbers in general case when the state is mixed and by 6 real 
numbers in case of the pure state. In the experiment we have to provide one more measurement in order to find the 
normalization. Working with single mode biphoton polarization states we can completely describe them by means 
of the Since the coherence matrix carry full information about the state it can be used as well as the density 
matrix. But density matrix calculations are more usual, so we are going to use them. Polarization density matrix can 
be found from correlation moments by the formula |18l llflj 
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/ 2A 2E 
V2D* C y/2F 
y 2E* V2F* 2B 



(3) 



Thus, after 9 measurements with the different orientations of the plates in the setup we have ehough information 
to reconstruct the density matrix of the polarization state. Moreover, this configuration of the measurement setup 
(fig. |3J) allows us to verify the orthogonality of the states that belong to the same basis. 



B. Compensation 

In order to have the three terms in superposition interfering, one must achieve their perfect overlap in frequency, 
momentum and time domains. From the experimental point of view this means that the biphoton wavepackets coming 
from the two type I crystals and from the type II crystal must be overlapped. The overlap in the frequency domain 
is achieved by the usage of 5nm bandwidth interference filters and the overlap in momentum is ensured by using 
pinholes that select one spatial mode of the biphoton field. But the overlap in time cannot be achieved easily when 
using a pulsed laser source, because it is necessary to compensate for all the group delays that biphoton wavepackets 
acquire during their propagation through the optical elements of the setup [20j . It was found that in order to overlap 
type-I biphotons with type-II, the pump pulse from the lower arm must be delayed. In our case the value of the delay 
is 50 ps. This was achieved by inserting an additional 2 mm BBO crystal in the lower arm. The overlap between 
the states |2, 0) and |0, 2) was achieved by inserting a 10 mm quartz plate directly after the two type I BBO crystals. 
After overlapping the biphotons with these techniques, the average coincidence count rate that we observed was of 
about 1 Hz. The high visibility of interference patterns that we obtained was a criterion for a good compensation (see 
figures Fig. Fig. [3 Fig. in next section). 



C. Calibration 



In order to create a given qutrit state we needed to have independent control over four real parameters - two relative 
amplitudes and two relative phases. 

Real amplitudes. In the experiment we used HWP1 to control the amplitude of the state |1,1), and HWP2 
to control the relative amplitudes of the states |2,0) and |0, 2). The calibration of these elements can be done by 
measuring moments related to pn, /O22, P33, defined in eq. (jjj in the tomography setup [Tt| . For our concrete 
experiment we had equal weights of each basis state in the superposition, so we had to set our waveplates in such 
a way that these moments fulfilled the condition: pn = P22 = P33- The values of the diagonal density matrix 
components were p u = 0.355 ± 0.032, p 22 = 0.34 ± 0.046, p 33 = 0.305 ± 0.035. 

Phase 0i3. The relative phase </>i3 = $3 — 4>i between states |2,0) and |0,2) can be controlled with the help 
of rotating quartz plates (QP2). The resulting state of biphoton field after these plates can be written as l^') = 
|2, 0) + e 2 ^ 13 10, 2). Varying the phase </>i3 by rotating the quartz plates (QP), one can observe an interference pattern 
when the polarizers in Brown-Twiss scheme are set to transmit +45° and —45° polarized light 15] . This effect, 
known as "space-time interference", can be used to calibrate $13, since we can assign value "0" to the position of the 
minimum of interference pattern, and "7r" to the position of the maximum. An obtained typical dependence from 
this parameter is presented at (Fig. 0J. A nonperiodic behavior is caused by a nonlinear dependence of an induced 
phase from the rotation angle. One can recalculate the dependence from the angle to dependence from the phase and 
the result is presented at (Fig. 0). This plot serves as a calibration curve for phase </>i3. 

It is important to note that in order to ensure the resulting state being close to a pure state it is necessary to 
achieve as high visibility as possible in space-time interference experiment with type I biphotons and polarization 
interference of type II biphotons (Fig. EJ) . The obtained visibilities of approximately 95% were considered as a good 
result in order to proceed with a final stage of the experiment - the combination of states |\E ,/ ) = |2,0) + e^ 13 |0, 2) 
and |1, 1) with a certain shift (f>\2 between them. 

Phase 4*12'. pump interference. The relation of the phase 4>i2 = 4>% — 4>x between the state |\&') and |1, 1) to the 
voltage applied on PZT can be found by monitoring the pump interference pattern in M-Z interferometer. We found 
that the change of voltage by 1 V results in phase shift of 0.33rad . 

Phase <fii2: generation of superposition between |2,0) and |1,1) states. We also measured an interference 
pattern for a state |\E'") = -^(|2,0) + exp(i0 12 )|l, 1)) varying the phase ^12. The polarization filters in Brown-Twiss 
scheme were set to measure moment A of the tomography protocol. Moment A exhibits no dependence on the phase 
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FIG. 4: Space-time interference pattern for state |^') = -^=(|2, 0) + e^ 13 j0,2)). Minimum corresponds to ^13 = it, maximum 
corresponds to cj>i3 = 0. 
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FIG. 5: Space-time interference pattern for state = t^(|2, 0) + e^ 13 10, 2}). Recalculated dependence on the phase (/>i2 
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FIG. 6: Polarization interference pattern from compensated type II biphotons. 



(j>i2 for a given state. Therefore we inserted an additional polarization transformer - a zero-order quarter waveplate 
operated at 800 nm and oriented at 20° with respect to the vertical direction. The unitary action of this waveplate 
is described by 3x3 matrix that converts an input state l^") to another state that is sent to a measurement scheme 
[T^ . For this certain waveplate, the state obtained after a transformer becomes 
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FIG. 7: Pump interference pattern in Mach-Zehnder interferometer. The distance between two maxima (in Volts) corresponds 

tO 012 = 7T. 



^ (-0.25 + 0.32i)exp(icf)i2) - 0.15 

0.41iearp(#i 2 ) + 0.25 + 0.32i | (4) 
^(0.25 + O.32i)exp(i0i 2 ) + 0.15 - 0.54z, 



and now moment A = c\c\ shows a certain dependence on 0i2. The experimentally obtained dependence (Fig.[5J) can 
also be used to calibrate the phase. It is well known that a period of this interference pattern coincides with the the 
pump interference period. The theoretical visibility of the interference is equal to 64%, while the visibility obtained 
in the experiment is 53%. We attribute this difference to the misalignments of the setup. 

We want to point out that in this experiment (for calibration purposes) we prepared a very specific superposition 
of basic states, interfering SPDC from type I and type II crystals, that was never reported in literature up to our 
knowledge. 




FIG. 8: The interference pattern for state l^ "*). The distance between two maxima (in Volts) also corresponds to 0i2 = it. 



D. Experimental procedure 

In order to create a given qutrit state we needed to have independent control over four real parameters - two relative 
amplitudes and two relative phases. In the experiment we used HWP1 to control the amplitude of the state |1, 1), 
and HWP2 to control the relative amplitudes of the states |2, 0) and |0, 2). The relative phase 0i3 = 03 — <pi between 
the states |2,0) and |0, 2) can be controlled with the help of rotating quartz plates (QP2). The relation of the phase 
4>i2 =02-01 between the state |*') = |2, 0) + e #13 |0, 2) and |1, 1) to the voltage applied to PZT can be found by 
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monitoring the pump interference pattern in M-Z interferometer. We found that the change of voltage by 1 V resulted 
in the phase shift of 51.7° and 4>i2 grew linearly with the applied voltage. 

States that constitute the first basis are trivial (Table I). They can be produced with the help of a single crystal, 
corresponding to type I or type II interaction. State |2, 0) is generated when first A/2 (HWP1) angle corresponds 
to the maximal reflection of the pump beam into the lower arm of a Mach-Zehnder and the angle of the second 
half-lambda waveplate (HWP2) is equal to 0°. In order to generate state |0,2), the HWP2 must be rotated by 45° 
degrees from 0°, and to generate state |1, 1) the HWP1 is rotated such, that the whole pump goes into the upper 
arm of Mach-Zehnder. Therefore, in the following, we will consider only the generation of the rest nine states, i.e. 
those forming the other three bases. According to Table I, only the relative phases between the basic states are to be 
varied. This allows us to use the same settings of the HWP's for the generation of nine states. It is also convenient to 
perform three sets of data acquisition - for the fixed 4>i3 values of 0, +120° and —120°, we change </>i2 values in the 
range of, say, few periods and perform all tomographic measurements for each value of the phase (f>i2 . Then we select 
the values of 4>\2 that correspond to the generation of the required state. For example, in order to generate the state 
/?', we use 0i3 = —120° and 0i2 = 120°. The values of the moments at this point allow us to restore a raw density 
matrix of the generated state and compare it to the theoretical value. 

The following procedure was used in order to verify the orthogonality of the states that form a certain basis. First 
we chose a set state to which we would tune our polarization filters. Then the values of the angles of quarter- and 
half- waveplates (Fig. |3J) (xi> X2, #2) that assure the maximal projection of the polarization state of each photon 
on the V direction can be calculated by mapping the set state on the Poincare sphere. Here, the lower index " 1" 
corresponds to the transmitted arm, and the index "2" to the reflected arm of BS We chose states \a'), \a") and \a'") 
to be our set states for each basis. Then, by setting the phase 13 fixed and by varying the phase 12 we measured 
the number of coincidence counts that correspond to the certain fourth order moment of the field. According to the 
orthogonality criterion, the coincidence rate should fall to zero when the values of 4>i% and 4>\2 correspond to the 
generation of the states orthogonal to the set ones. 



E. Results and discussion 



Let us consider the generation of the state |/3")- In this case ^13 = 0, $12 = 120. In Fig. [5]the measured values of 
the real and imaginary parts of the density matrix components P21 and ^32 on phase </>i2 are shown as function of the 
phase 012- The number of accidental coincidences was negligibly small and was not subtracted in data processing. 
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FIG. 9: Imaginary and real values of non-diagonal density matrix components used to reconstruct state \f3"). Theoretical 
dependence is plotted with a solid curve. 



The phase 0i 3 = remained constant during the tomography procedure. After obtaining the dependence of 
the moments P21 and P32 on phase (f>i2 we fitted our data with theoretical dependencies, using the least-square 
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approximation method. The obtained values of all components were substituted in Eq. The obtained density 
matrix for state \/3") is given below. 



P/3" 



/ 0.355 -0.054 - 0.210J 0.315 - O.OlOi 

-0.054 + 0.210i 0.340 -0.106 + 0.262j 

\ 0.315 + 0.010? -0.106 - 0.262i 0.305 



(5) 



The eigenvalues of this matrix are Ai = 0.877, A2 = 0.136, A3 = —0.013. A corresponding set of eigenvectors 
is X = (0.587,-0.173 + 0.52K, 0.594 - 0.07K);Y = (0.642, 0.379 - 0.649i, 0.048 + 0.143z); Z = (0.493,-0.287 + 
0.224z, —0.769 — 0.178z). Although the density matrix (Eq. [SJ is Hermitian and the condition Tr(p) = 1 is satisfied, 
it doesn't correspond to any physical state because of the negativity of one of the eigenvalues. We want to point out 
that a first main component (p\ xp )ij = XiX* of a considered density matrix, which has a weight 0.878 is already 
close to the theoretical state vector \[3") and the corresponding fidelity is F = Tr{p t hP 1 exp ) = 0.9903. The other 
two components correspond to the "experimental noise" that is due mainly to misalignments of a setup and small 
volume of collected data. We have obtained similar eigenvalues for all other states and raw fidelity computed for the 
main density matrix component as described above have varied from 0.983 to 0.998. We also employed the maximum 
likelihood method of quantum state root estimation (MLE) 0, to make a tomographically reconstructed matrix 
satisfy its physical properties, such as positivity. The results are presented in the following table (Table II). The 
level of statistical fluctuations in fidelity estimation was determined by the finite size of registered events (~ 500). 
All experimental fidelity values lie within the theoretical range of 5%(F = 0.9842) and 95%(F = 0.9991) quantiles 
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TABLE III: Fidelities estimated with Maximum Likelihood Method. 



The obtained fidelity values show the high quality of the prepared states. The other test of the quality of states 
is the fulfillment of the orthogonality criterion for the states that belong to the same basis. For each set state we 
calculated the settings of waveplates in our measurement setup that ensured the maximal projection of each photon 
on the vertical polarization direction. In Fig. we show the dependence of the coincidence rate for the following 
setting of waveplates %i = 28.3°, X = -33.5°, \2 = -24°,6> 2 = -2°. 




FIG. 10: Dependence of number of coincidences on a phase <f>\2 for a given settings of polarization filters. 



These values correspond to the set state \a"'). As one can see, for the fixed value 4>i3 = the coincidence rate 
is almost equal to zero, when phase 0i2 = —120°. This corresponds to the generation of the state \(3"'), which is 
orthogonal to \a"'). The visibility of this pattern is equal to 93.2%. For the other bases, the obtained values of 
visibilities varied from 92% to 95%. With these values of visibility, the lowest value of coincidence rate corresponds 
to the accidental (Poissonian) coincidence level and therefore the obtained data verifies the orthogonality criterion. 
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Here we try to clarify how we understand the role of statistics in quantum state reconstruction. When we find that 
fidelity F = 0.990, then this result is sufficient without pointing out +/-0.XXX, because the obtained value shows 
only the degree of correspondence between the desirable result (to obtain fidelity as close to unity as possible) and 
the achieved precision of quantum state reconstruction. As usual, all works on quantum state tomography end up on 
this. At the same time we think that if wc consider the question more thoroughly, it is necessary not only to point out 
the result, but also to try to answer the question about the principal precision that we can obtain in the experiment 
if we consider the finite (and not so large) volume of statistical data. Luckily we can also answer this question; it is 
overviewed in the paper ^| . On the basis of theory of statistical fluctuations of state vector estimation, that was 
developed in , we obtained that the interval of the expected statistical fidelity fluctuations lies within 0.9842 (5% 
quantile) and 0.9991 (95% quantile). 

The fact that all obtained fidelity values lie within that interval shows the correspondence of the experiment to the 
statistical theory. We also point out that due to the complexity of experiment, the volume of the data is not large 
(about 500 events in each of nine experiments). If we increase the volume of data, then, sooner or later, it will come to 
the point, when the correspondence of experiment to statistical theory will be violated. The experimentally achieved 
fidelity value would be lower then the theoretical due to the inevitable instrumental errors of experiment (see also 
[l8j ') . In this case the quality of experiment can be estimated by the so called "coherence volume" of data. If we 
exceed that volume, then the fidelity "saturates" on the level which is somewhat less then unity. All further increase 
of the volume will not lead to the increase of fidelity since the precision of state reconstruction will be limited by the 
instrumental errors and instability of the experimental setup. 

III. CONCLUSIONS 

We realized an interferometric method of preparing the three-level quantum optical systems, that relied on the 
polarization properties of single-mode two-photon light. The specific sequence of states was generated and measured 
with high fidelity values. The orthogonality of the states that form mutually unbiased bases was experimentally 
verified. As an advantage of this method we note that all control of the amplitudes and phases of each basic state 
in superposition Q is done using linear optical elements, making it easy to switch from one state to another and 
providing the full control over the state (JU . The main disadvantage is that we cannot generate an entangled qutrits 
in this configuration. Our setup also allows one to prepare an arbitrary polarization qutrit state on demand. 
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